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1 .  INTRODUCTION 


1.1  Background .  The  work  being  reported  is  part  of  a  current  research 
task  at  the  Ballistic  Research  Laboratory  (BRL)  to  design  a  large  blast/ 
thermal  simulator  for  multiservice  use.  The  simulator  is  to  be  large 
enough  to  conduct  full-size  vehicle  tests  throughout  a  broad  spectrum  of 
blast/thermal  effects,  simulating  a  range  of  blast  and  thermal  from  nuclear 
environments.  The  simulator  facility  is  to  be  similar  to  the  French  facil¬ 
ity  of  Gramat  (Reference  1),  but  is  to  have  expanded  operational  capabili¬ 
ties,  both  in  size  and  in  the  simulation  ranges.  Previous  phases  of  the 
task  at  BRL  have  been  reported  in  Reference  2,  3,  and  A. 

1.2  Objectives .  The  experiments  reported  here  were  designed  to  answer 
some  questions  concerning  the  basic  shock  parameters  necessary  to  the 
design  of  the  full-size  blast  simulator.  These  parameters  are:  the  test 
station  pressure  as  a  function  of  driver  pressure,  the  driver  length/volume 
needed  to  produce  the  desired  range  of  yields,  the  cold  gas  driver  effect, 
and  the  effect  of  the  recompression  fan  from  the  expansion  of  the  divergent 
throat  nozzle.  These  effects  were  to  be  examined  both  experimentally  and 
by  hydrocode  simulation.  References  2,  A,  and  5  describe  the  BRL-Q1D 
hydrocode  used  for  the  computer  simulation. 

2.  TEST  PROCEDURES 

The  experimental  parameter  study  task  consisted  first,  of  the  design 
and  fabrication  of  a  1:57  scale,  single  cold  air  driver  model  of  the  antic¬ 
ipated  full-size  blast  simulator;  secondly,  of  a  selection  of  a  range  of 
driver/throat  baffle  configurations  to  produce  the  desired  test  pressure/ 
yield  levels,  and  to  record  the  pressure-time  results  as  a  function  of  test 
section  location.  The  model  is  described  in  the  next  section. 

2.1  Simulator  Model.  The  1:57  scale  for  the  model  was  chosen  for  the 
experiments,  primarily  because  of  the  ready  availability  of  the  steel  pipe 
components  of  the  simulator  model.  A  sketch  of  the  shock  tube  model  is 
shown  in  Figure  1.  The  drivers  consisted  of  assorted  lengths  (6.05-288.0 
cm)  of  10.16  cm  inside  diameters,  smooth  thick-walled  pipe.  This  was 
attached  to  the  long  test  section  (.25A  I.D.  by  17. 1A  m  long  -  closed  at 
end)  through  a  converging  throat  or  baffle  to  the  diaphragm  section.  A 
diverging  section  expanded  to  near  the  inside  diameter  of  the  test  section. 
Cold  compressed  air  was  used  for  the  driver  gas.  Mylar,  aluminum,  and 
copper  diaphragms  were  used  to  contain  the  monitored  driver  pressure  until 
the  diaphragm  self-ruptured . 

Stations  for  pressure  transducers  (side-on  and  stagnation  probes)  were 
located  along  the  test  section  at  7,  15,  20,  2A,  and  28  test  section 
diameters  as  noted  on  the  drawing  in  Figure  1.  The  stations  are  listed  on 
the  pressure-time  plots  as  70-280,  respectively,  in  the  Results  section. 

The  recording  and  data  reduction  instrumentation  is  described  in  the 
next  section. 

2.2  Instrumentation .  Two  separate  sets  of  recording  instrumentation 
were  used  for  the  experimental  shots.  Figures  2-A  and  2-B  describe  the  two 
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NOTES: 

(1)  DIMENSIONS  IN  CENTIMETERS 

(2)  NOT  TO  SCALE 

(3)  DIAPHRAGMS  WERE  MYLAR,  AL  or  Cu 

(4)  DIMENSIONS  ARE  FOR  INSIDE 


Figure  1.  Sketch  of  1:57  Scale  Single  Driver  Model  of  Large  Blast  Simulator 
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A.  Tape  System 


Figure  2.  Schematic  of  Data  Acquisition/Reduction  Systems. 
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B.  Oscillograph  System 


Figure  2.  Schematic  of  Data  Acquisition/Reduction  Systems  (Cont'd). 


4 


types  used.  Both  were  used  to  record  the  output  generated  by  the  PCB  Model 
113A24  quartz  pressure  transducers  (Reference  6).  The  tape  recording 
systems  (Figure  2-A)  was  used  for  the  initial  multi-station  recording. 

Later  in  the  test  program  only  two  stations  were  used,  so  the  second  system 
(Figure  2-B)  with  the  digitizing  oscilloscope  was  used. 

Both  systems  recorded  the  output  from  the  transducers  coupled  through 
power  supply  interface  cards  and  data  amplifiers  to  the  recorders.  On-site 
comparisons  of  the  pressure-time  plots  were  made  directly  from  the  hard 
copies  produced.  Final  data  processing  was  completed  with  the  computer, 
printer,  and  plotter.  Plots  of  pressure-time  records  for  various  driver 
configurations  are  shown  in  the  Results  section. 

2.3  Test  Matrix.  A  variety  of  shots  were  fired  with  the  various 
driver  configurations.  The  test  matrix  is  given  in  Table  1  for  several 
shot  conditions.  The  shots  are  listed  by  number,  with  driver  parameters, 
ambient  test  conditions,  and  test  section  overpressure  at  Station  70  (7 
dia.)  listed  for  reference. 

The  driver  lengths  varied  from  a  minimum  of  6.05  cm  to  288  cm  for  the 
maximum  "long"  driver.  The  inside  diameter  of  the  drivers  remained  con¬ 
stant  at  10.16  cm.  The  pressure  in  the  driver  varied  from  about  300  kPa  to 
a  little  over  18000  kPa.  Throat  baffles  at  the  diaphragm  section  varied 
from  a  test  section  to  throat  baffle  area  ratio  of  16.1  to  64:1.  Sam¬ 
ple  pressure-time  records  and  a  discussion  of  their  characteristics  are 
given  in  the  Results  section. 

3.  RESULTS 

The  primary  results  from  the  test  program  are  the  pressure-time  records 
and  their  impulses  from  the  test  stations  for  the  various  test  configura¬ 
tions.  This  section  will  discuss  the  change  of  pressure-time  records 
(side-on,  stagnation, and  dynamic)  as  the  driver  parameters  were  changed. 
Stations  at  7  and  20  diameters,  in  particular,  will  be  compared  with  com¬ 
puter  code  predictions.  Both  short  (11.16  cm)  and  long  (288  cm)  drivers 
will  be  used  for  the  comparisons.  Finally,  particular  examples  of  cold 
gas/recompression  fan  effects  will  be  shown. 

3.1  Test  Station  Pressure  as  a  Function  of  Driver  Parameters.  The 
first  set  of  records  shown  in  Figure  3-6  illustrate  the  decay  of  peak  pres¬ 
sure  along  the  test  section  for  a  short  driver  (11.12  cm)  as  a  function  of 
driver  pressure.  The  driver  pressure  was  varied  within  the  range  314  kPa 
to  14479  kPa  for  this  series  of  shots.  These  values  are  plotted  in  the 
graph  of  Figure  7.  Initial  test  section  pressures  (20-225  IcPa)  measured  at 
7  diameters  decay  with  distance  along  the  test  section  to  a  station  at  28 
diameters.  Decays  of  about  25  to  43%  were  measured  over  this  distance. 

The  hydrocode  predicted  similar  decays,  but  for  driver  pressure  above 
3000-4000  kPa  the  predicted  shock  overpressure  was  higher  than  the  experi¬ 
mental  values .  Below  this  driver  range  the  predicted  values  varied  both 
above  and  below  the  experimental  values . 

It  can  be  seen  from  the  family  of  curves  plotted  from  the  long  driver 
data  and  shown  in  Figure  8,  that  there  is  very  little  attenuation  of  the 
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The  text  is  continued  on  page  21 . 


TABLE  1.  TEST  MATRIX 


Shock 


Shot 

Driver 

Driver 

Test  Section  to 

Ambient 

Ambient 

Overpressure 

Number 

Pressure 

Length 

Throat  Area  Ratio 

Pressure 

Temperature 

Station  70 

kPa 

cm 

kPa 

°k 

kPa 

4 

314 

11.13 

16:1 

102.4 

294.8 

20.0 

5 

617 

102.3 

295.0 

32.6 

6 

1124 

102.2 

295.2 

49.3 

7 

2055 

102.1 

296.6 

75.0 

S 

3137 

102. 1 

297.1 

96.6 

9 

4413 

Same 

Same 

101.9 

296.0 

118.2 

10 

5240 

103.2 

292.9 

137.0 

11 

7757 

103.3 

294.2 

169.6 

12 

11273 

103.3 

294.2 

188.0 

13 

14479 

103.3 

295.4 

225.0 

14 

5171 

288 

16:1 

102.9 

295.2 

131.5 

15 

7481 

102.2 

297.9 

155.0 

16 

3103 

102.2 

297.7 

88.0 

17 

603 

102.2 

297.7 

23.5 

18 

300 

102.2 

298.1 

13.3 

19 

1138 

Same 

Same 

102.0 

297.7 

38.0 

20 

2034 

101.9 

297.5 

60.4 

21 

4220 

102.0 

297.0 

106.5 

22 

10963 

102.0 

297.3 

193.0 

23 

13169 

101.7 

297.7 

215.0 

29 

1103 

11.13 

64:1 

102.5 

297.2 

29.2 

30 

2034 

102.5 

297.4 

39.1 

31 

3103 

102.3 

298.2 

56.6 

32 

8618 

102.5 

297.5 

107.0 

TABLE  1.  TEST  MATRIX  (CONT'D) 


Shot 

Driver 

Driver 

Test  Section  to 

Ambient 

Ambient 

Shock 

Number 

Pressure 

Length 

Throat  Area  Ratio 

Pressure 

Temperature 

Overpressure 

kPa 

cm 

kPa 

°k 

kPa 

33 

565 

11.13 

33:1 

102.2 

297.9 

26.1 

34 

1069 

102.2 

298.3 

36.2 

35 

2027 

102.2 

297.5 

54.3 

36 

3130 

102.2 

297.6 

74.2 

37 

5192 

Same 

Same 

102.2 

298.1 

99.6 

38 

8618 

102.2 

298.8 

137.2 

39 

10515 

102.6 

298.9 

156.2 

40 

14479 

102.6 

298.6 

195.3 

45 

13617 

26.36 

16:1 

102.1 

298.6 

217.9 

46 

13962 

33.98 

16:1 

102.0 

298.9 

205.5 

47 

317 

101.1 

297.5 

18.3 

48 

579 

101.1 

297.2 

27.0 

49 

1096 

101.0 

296.5 

43.5 

50 

2054 

101.0 

296.5 

76.0 

51 

607 

67.0 

64:1 

102.5 

298.5 

16.3 

52 

690 

89.86 

64:1 

102.5 

298.5 

22.8 

53 

338 

16:1 

101.7 

297.6 

14.6 

54 

710 

64:1 

101.5 

297.1 

15.5 

55 

359 

178.76 

16:1 

101.4 

297.0 

14.3 

56 

359 

206.70 

102.2 

296.7 

13.3 

57 

352 

237.18 

102.2 

297.5 

13.2 

58 

614 

125.42 

64:1 

102.2 

296.6 

10.0 

59 

324 

16:1 

102.2 

297.9 

17.9 

60 

703 

64:1 

102.6 

296.7 

10.0 

61 

896 

33:1 

102.6 

296.7 

20.0 

62 

717 

33:1 

102.5 

296.7 

20.0 

63 

510 

33:1 

102.7 

296.7 

15.0 

TABLE  Is  TEST  MATRIX  (CONT'D) 


Shot 

Driver 

Driver 

Test  Section  to 

Ambient 

Ambient 

Shock 

Number 

Pressure 

Length 

Throat  Area  Ratio 

Pressure 

Temperature 

Overpressure 

kPa 

cm 

kPa 

°k 

kPa 

64 

17754 

33.98 

16:1 

102.6 

296.4 

262.0 

65 

17582 

6.05 

102.5 

296.4 

200.0 

66 

13962 

18.74 

102.2 

298.0 

222.7 

67 

13962 

26.36 

Same 

102.1 

297.9 

224.8 

68 

14307 

33.98 

102.0 

298.0 

237.5 

69 

3137 

102.1 

298.0 

106.5 

70 

5171 

102.1 

298.0 

147.2 

71 

8446 

102.2 

298.1 

180.0 

72 

11549 

102.2 

298.0 

225.0 

73 

13445 

102.57 

103.0 

298.0 

220.0 

74 

13169 

72.09 

103.0 

298.4 

221.3 

75 

3137 

72.09 

103.0 

297.9 

110.0 

77 

3123 

110.18 

16:1 

103.0 

298.1 

106.0 

78 

1117 

110.18 

102.9 

297.8 

47.2 

79 

1124 

72.09 

102.7 

297.9 

48.5 

80 

1145 

145.74 

102.6 

297.6 

49.5 

81 

483 

33.99 

33:1 

101.8 

298.4 

19.3 

32 

1827 

101.7 

298.5 

47.5 

83 

5199 

101.7 

298.4 

103.0 

84 

14789 

101.7 

298.4 

200.2 

85 

483 

67.01 

Same 

101.6 

298.2 

20.0 

86 

1827 

101.6 

297.9 

46.8 

87 

5240 

101.6 

297.9 

108.0 

83 

14651 

101.5 

297.9 

200.0 

89 

496 

94.95 

101.7 

298.7 

14.6 

90 

1834 

94.95 

101.6 

297.9 

40.8 

91 

545 

145.75 

101.6 

298.0 

15.0 

92 

490 

11.13 

101.8 

297.9 

19.5 

93 

18271 

6.05 

101.9 

298.4 

172.5 

94 

903 

74.62 

64:1 

102.7 

298.0 

13.8 

KPA 


f<PA-KSEC  IMPULSE.  KPA-NSEC  IMPULSE.  KPA-MSEC  IMPULSE.  KPA-MSEC 


TEST:  LBTS/57 


TEST:  LBTS/57 


TEST;  LBTS/57 


TEST:  LBTS/57 


B.  Stagnation  Overpressure 


Figure  4. 


Pressure-Time  Records  from  Test  Section 


Short  Driver,  1124  kPa  (Cont'dl. 
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IMPULSE.  KPA-HSEC  IMPULSE.  KPA-MSEC  IMPULSE.  KPA-MSEC  IMPULSE.  KPA-NSEC 


PRESSURE,  KPA 


TEST.  LBTS/F7 


0  5  10  15  20 


MSEC 


TEST.  LBTS/57 


MSEC 


A.  Side-On  Overpressure 


Figure  5.  Pressure-Time  Records  from  Test  Section  -  Short  Driver,  4413  kPa. 


KPA-MSEC  IMPULSE.  KPA-MSEC  IMPOSE.  KPA-MSEC 


PRESSURE,  KPA 


TEST:  LBTS/57 


0  5  10  15  20 

TDC,  MSEC 


TEST:  LBTS/57 


TDE,  MSEC 


TEST:  LBTS/57 


0  5  10  15  20 


TIME,  MSEC 


TEST:  LBTS/57 


TOE,  MSEC 

A.  Side-On  Overpressure 


Pressure-Time  Records  from  Test  Section 


Short  Driver,  14479  kPa. 
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KPA-HSEC  IMPULSE.  KPA-HSEC  IMPULSE.  KPA-HSEC  IMPULSE.  KPA-MSEC 


TEST:  LBTS/57 


0  5  10  15  20 

TD€,  MSEC 


TEST-  LBTS/57 


TDfL  MSEC 


TEST:  LBTS/57 


0  5  10  15  20 


TDC,  MSEC 


TEST.  LBTS/57 


TIME,  MSEC 


B.  Stagnation  Overpressure 


Figure  6.  Pressure-Time  Records  from  Test  Section  -  Short  Driver,  14479  kPa  (Cont'd). 
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.  KPA— MSEC  IMPULSE.  KPA-HSEC  Itf’ULSE.  KPA-KSEC  imJLSE.  KPA-MSEC 


SHOCK  OVERPRESSURE,  kPa 


Figure 


SINGLE  DRIVER  -  10.16cm  DlA.x  288cm  LONG, 


8.  Side-On  Shock  Overpressure  Along  Test  Section  -  Long  Drivi 


2 


overpressure  in  traveling  from  7  to  28  test  section  diameters  since  "flat- 
top"  shock  waves  have  only  viscous  attenuation.  The  hydrocode  has  no 
program  for  viscous  attenuation  and,  therefore,  shows  no  decay  at  all. 

Figure  9  is  an  operating  curve  showing  the  peak  overpressure  to  be 
expected  at  the  two  important  stations  (7  and  20  dia.)  as  a  function  of  the 
driver  pressure.  The  values  plotted  in  Figures  8  and  9  were  taken  from  the 
pressure-time  records  shown  in  Figures  10-13.  Typical  examples  of 
pressure-time  records  of  side-on  overpessure,  stagnation  overpressure,  and 
their  subtraction  (the  dynamic  pressure)  are  shown  for  the  long  driver  of 
288  cm. 


Generally,  between  3000  and  4000  kPa  driver  pressure,  the  predicted 
(computer  code  simulation)  and  experimental  shock  overpressures  deviated. 

At  the  extreme  high  end  for  13500  kPa,  the  predicted  value  of  pressure  was 
about  270  kPa.  Experimentally,  the  value  was  215  kPa .  The  lower  experi¬ 
mental  values  may  be  explained  by  the  thicker  diaphragm  needed  at  the 
higher  driver  pressures.  Opening  times  are  finite.  Obviously,  there  is 
not  a  complete  diaphragm  removal  at  zero-time  as  was  assumed  by  the  hydro¬ 
code.  Incomplete  breaking  of  the  diaphragm  also  occurred.  All  these  fac¬ 
tors  tended  to  produce  lower  test  section  pressures  than  predicted  from  the 
BRL-Q1D  code  for  the  3000-4000  kPa  range.  Below  this  range,  the  predicted 
values  varied  above  and  below  the  experimental  values . 

Tables  2-4  summarize  the  experimental  data  according  to  driver  length 
and  pressure.  Representative  values  of  the  blast  parameters  are  given  for 
the  two  test  stations  at  7  and  20  diameters.  Additionally,  scaled-up 
values  of  the  parameters  for  the  full-size  simulator  are  given  in  the  last 
four  columns  of  the  tables.  These  values  were  obtained  by  multiplying  by 
the  scale  factor  of  57.  The  last  column  of  the  tables  display  the  equiva¬ 
lent  high  explosive  yield  for  the  full-size  conditions.  The  scaling  cal¬ 
culations  will  be  discussed  in  some  detail  in  the  Analysis  section. 

3.2  Hydrocode  Results.  One  of  the  basic  research  efforts  at  BRL  is  to 
computationally  model  shock  tube  processes.  The  BRL-Q1D  code  used  for 
these  predictions  is  an  adiabatic,  inviscid  eulerian  computer  algorithm 
adapted  at  BRL  for  this  purpose.  See  References  2,  4,  5,  and  7  for  addi¬ 
tional  descriptions  and  uses  of  the  code.  All  references  to  hydrocode  com¬ 
putations  refers  to  the  BRL  Quasi-one-dimensional  code  (BRL-Q1D) . 

The  one-dimensional  computational  modeling  used  had  the  advantages  that 
it  was  inexpensive  to  run,  had  capability  to  perform  parametric  stud¬ 
ies  quickly,  and  was  useful  for  obtaining  good  engineering  approximations 
to  the  simulator  design.  Its  obvious  disadvantage  was  that  the  one¬ 
dimensional  analysis  was  necessarily  an  approximation.  As  is  seen  later, 
the  BRL-Q1D  code  did  predict  appearance  of  the  cold  driver  gas  and  the 
recompression  fan  from  the  diverging  nozzle,  but  mispredicted  the  magnitude 
of  both  effects  and  their  arrival  after  the  shock  front  at  a  particular 
test  station. 

Figures  14-17  show  typical  records  obtained  for  the  short  driver  (11.13 
cm).  Stations  1,  2,  and  3  of  the  computer  runs  correspond  to  Stations  70, 
150,  and  200  in  the  experimental  records.  As  noted  above  in  Section  3.1  of 
the  Results,  the  code  predicted  higher  peak  overpressure  in  the  test  sec¬ 
tion  for  driver  pressures  in  the  3000-4000  kPa  range.  Below  this  range, 


2i  The  text  is  continued,  on  page  43. 


DRIVER  OVERPRESSURE,  k 


Figure  9.  Test  Station  Overpressure  as  a  Function  of  Driver  Pressure 


PRESSURE,  KPA 


g  150 


nti! 


STATION-  158 


STATION:  200 


IMPULSE.  KPA-MSEC  IMPULSE.  KPA-HSEC  M'LLSE.  KPA-MSEC  DRLSI.  KPA-MSEX 


PRESSURE,  KPA 


IMPULSE.  KPA-MSEC  IMPULSE.  KPA-NSEC  IhPULSE.  KPA-MSEC  IMPULSE.  KPA-MSEC 


TEST:  LBTS/57 
SHOT.  15 

STATION:  150  8TA8 


PRESSURE,  KPA 


IMPULSE.  KPA-MSEC  IMPULSE.  KPA-MSEC  BPULSE.  KPA-MSEC  IMPULSE.  KPA-MSEC 


TEST:  LBTS/57 


TEST:  LBTS/57 


TEST:  LBTS/57 


TEST:  LBTS/57 


Figure  13.  Pressure -T i me  Records  from  Test  Section  -  Long  Driver,  15169  kPa  (Cont'd). 
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IMPULSE.  KPA-KSEC  IMPULSE.  KPA-MSEC  IMPULSE.  KPA-MSEC  IMPULSE.  KPA-MSEC 


TABLE  2.  BLAST  WAVE  PARAMETERS  VS  DRIVER  PARAMETERS 
TEST  SECTION  TO  THROAT  AREA  RATIO  OF  16:1 


Shot 

Driver 

Driver 

Side-On 

Number- 

Pressure 

Length 

Overpressure 

Station 

kPa 

cm 

kPa 

4-70 

314 

11.13 

20.0 

200 

17.0 

5-70 

617 

33.1 

200 

24.9 

6-70 

1124 

49.9 

200 

same 

34.1 

8-70 

3137 

100.1 

200 

66.7 

13-70 

14479 

226.2 

200 

155.8 

67-70 

13962 

26.36 

224.8 

200 

209.6 

47-70 

317 

33.98 

18.6 

200 

16.5 

49-70 

1096 

43.4 

200 

49.1 

69-70 

3137 

106.2 

200 

94.6 

72-70 

11549 

224.8 

200 

202.5 

Full  Size 


Positive 

Positive 

Driver 

Driver 

Positive 

Duration 

Impulse 

Length 

Volume 

a 

Impulse 

ms 

kPa-ms 

m 

3 

m 

kPa-s 

3.7 

22.5 

6.34 

166.93 

1.28 

4.5 

22.8 

1.30 

4.7 

45.0 

2.57 

5.8 

45.0 

2.57 

6.0 

75.1 

4.28 

6.9 

77.3 

4.41 

5.1 

170.1 

9.70 

7.3 

185.5 

10.57 

7.6 

573.9 

32.71 

12.2 

710.9 

40.52 

8.5 

803.9 

15.03 

395.74 

45.82 

18.0 

1249.3 

71.21 

8.4 

54.5 

19.37 

510.01 

3.11 

8.8 

56.1 

3.20 

13.8 

173.5 

9.89 

16.5 

168.7 

9.62 

12.3 

376.7 

21.47 

15.0 

421.4 

24.02 

12.0 

830.4 

47.33 

20.3 

1317.5 

75.09 

2 

^Yields  are  based  on  full  size  driver  diameter  of  5.79  m,  area  of  26.33  m  ,  and  are  calculated 
from  side-on  overpressure  impulse.  Last  four  columns  are  for  full  size  for  Tables  2-4. 


Yield* 

kt 

0.08 

0.12 

0.27 

0.36 

0.65 

1.29 

3.0 

6.3 
67.0 

161.0 

185.0 

732.0 

1.4 

1.9 

9.9 
7.6 

31.0 

49.0 

205.0 

882.0 


TABLE 

2 .  BLAST  WAVE 

PARAMETERS  VS  DRIVER 

PARAMETERS ; 

TEST 

SECTION  TO 

THROAT  AREA 

RATIO  OF 

16:1  (C0NT.; 

) 

Full 

Size 

Shot 

Driver 

Driver 

Side-On 

Positive 

Positive 

Driver 

Driver 

Positive 

Number- 

Pressure 

Length 

Overpressure 

Duration 

Impulse 

Length 

Volume 

Impulse 

Yield 

Station 

kPa 

cm 

kPa 

ms 

kPa-ms 

m 

3 

m 

kPa-s 

kt 

79-70 

.  1124 

72.09 

48.5 

29.0 

351.7 

41.09 

1081.89 

20.04 

74.0 

200 

42.5 

30.5 

359.6 

20.50 

91.0 

75-70 

3137 

110.0 

— 

850.6 

48.48 

344.0 

200 

94.0 

— 

857.5 

48.88 

411.0 

74-70 

13169 

221.3 

9.5 

1319.5 

75.21 

831.0 

200 

233.0 

30.8 

2571.3 

146.56 

5869.0 

53-70 

338 

89.86 

13.5 

19.2 

118.0 

51.22 

1348.60 

6.73 

27.0 

200 

14.0 

20.2 

131.0 

7.47 

34.0 

73-70 

13445 

102.57 

220.0 

10.3 

1412.0 

58.46 

1539.24 

80.48 

1022.0 

200 

224.8 

43.2 

3800.0 

216.60 

19553.0 

78-70 

1117 

110.18 

47.2 

41.4 

512.2 

62.81 

1653.77 

29.20 

232.0 

200 

40.5 

44.8 

505.9 

28.84 

278.0 

77-70 

3123 

106.0 

— 

1350.0 

76.95 

1424.0 

200 

99.0 

57.3 

1274.8 

72.66 

1275.0 

80-70 

1145 

144.97 

49.5 

57.0 

700.0 

83.07 

2187.20 

38.18 

473.0 

200 

45.0 

56.0 

720.0 

41.04 

677.0 

59-70 

324 

125.42 

17.9 

26.0 

173.1 

71.49 

1882.32 

9.87 

49.0 

200 

16.6 

39.4 

180.7 

10.30 

64.0 

55-70 

359 

178.76 

14.2 

37.3 

265.0 

101.93 

2683.79 

15.11 

275.0 

200 

13.0 

38.5 

272.8 

15.55 

365.0 

57-70 

352 

237.18 

13.2 

48.7 

328.8 

135.19 

3559.52 

18.74 

618.0 

200 

13.0 

48.5 

327.4 

18.66 

632.0 

TABLE  3.  BLAST  WAVE  PARAMETERS  VS  DRIVER  PARAMETERS; 

TEST  SECTION  TO  THROAT  AREA  RATIO  OF  33:1 


Full  Size 


Shot 

Driver 

Driver 

Side -On 

Positive 

Positive 

Driver 

Driver 

Positive 

Number - 

Pressure 

Length 

Overpressure 

Duration 

Impulse 

Length 

Volume 

Impulse 

Yield 

Station 

kPa 

cm 

kPa 

ms 

kPa-ms 

m 

3 

m 

kPa-s 

kt 

93-70 

18271 

6.05 

172.5 

7.3 

523.1 

3.45 

90.84 

29.82 

61.0 

200 

135.0 

11.3 

558.3 

31.82 

84.0 

92-70 

490 

11.13 

19.5 

6.0 

30.5 

6.34 

166.93 

1.74 

0.22 

200 

15.0 

7.5 

32.1 

1.83 

0.44 

35-70 

2027 

54.3 

11.3 

120.0 

6.84 

2.3 

200 

39.0 

13.8 

120.0 

6.84 

4.0 

37-70 

5192 

99.6 

9.1 

232.5 

13.25 

7.7 

200 

83.4 

11.0 

271.0 

15.45 

14.9 

40-70 

14479 

195.3 

8.3 

567.9 

32.37 

72.0 

200 

142.4 

14.4 

636.5 

36.28 

121.0 

81-70 

483 

33.99 

19.3 

20.8 

88.0 

19.37 

510.01 

5.02 

5.5 

200 

15.7 

18.5 

84.4 

4.81 

7.3 

82-70 

1827 

47.5 

33.8 

284.9 

16.24 

39.0 

200 

39.3 

35.3 

293.7 

16.74 

58.0 

83-70 

5199 

103.0 

43.5 

686.2 

39.11 

192.0 

200 

88.9 

36.3 

659.6 

37.60 

198.0 

84-70 

14789 

200.2 

17.7 

937.2 

53.42 

321.0 

200 

200.0 

25.5 

1525.6 

89.96 

1386.0 

85-70 

483 

67.01 

20.0 

34.2 

153.1 

38.19 

1005.53 

8.73 

27.0 

200 

15.7 

34.3 

159.5 

9.09 

51.0 

86-70 

1827 

46.8 

80.5 

560.0 

31.92 

304.0 

200 

40.6 

63.3 

550.0 

31.35 

356.0 

87-70 

5240 

108.0 

— 

1275.0 

72.68 

1181.0 

200 

94.9 

— 

1300.0 

74.10 

1415.0 

88-70 

14651 

191.7 

15.0 

1450.0 

82.65 

1218.0 

200 

200.0 

41.0 

2640.0 

150.48 

7180.0 

89-70 

496 

94.95 

14.6 

47.0 

200.0 

54.12 

1424.97 

11.40 

114.0 

200 

13.5 

47.5 

202.0 

11.51 

138.0 

90-70 

1834 

40.8 

80.0 

726.0 

41.38 

804.0 

200 

39.8 

— 

720.0 

41.04 

771.0 

63-70 

510 

125.42 

15.0 

39.4 

277.9 

71.49 

1882.32 

15.84 

288.0 

200 

13.8 

67.3 

314.9 

17.95 

508.0 

91-70 

545 

145.75 

15.2 

69.0 

315.0 

83.07 

2187.21 

17.96 

414.0 

200 

15.0 

72.3 

363.0 

20.69 

641.0 

TABLE  4.  BLAST  WAVE  PARAMETERS  VS  DRIVER  PARAMETERS; 

TEST  SECTION  TO  THROAT  AREA  RATIO  OF  64:1 


CrJ 

oo 


Shot 

Driver 

Driver 

Side-On 

Positive 

Number- 

Pressure 

Length 

Overpressure 

Duration 

Station 

kPa 

cm 

kPa 

ms 

29-70 

1103 

11.13 

29.2 

28.5 

200 

21.8 

33.8 

30-70 

2034 

39.1 

21.6 

200 

28.8 

33.2 

31-70 

3103 

56.6 

25.8 

200 

39.6 

34.0 

32-70 

8618 

107.0 

42.1 

200 

86.9 

35.1 

94-70 

903 

74.62 

13.8 

>85.0 

200 

13.0 

>85.0 

54-70 

710 

89.86 

15.5 

86.0 

200 

13.8 

86.0 

Full  Size 


Positive 

Driver 

Driver 

Positive 

Impulse 

Length 

Volume 

Impulse 

Yield 

kPa-ms 

m 

3 

m 

kPa-s 

kt 

76.8 

6.34 

166.93 

4.38 

1.8 

84.0 

4.79 

4.1 

112.8 

6.53 

3.6 

134.0 

7.64 

9.7 

187.5 

10.69 

10.4 

177.0 

10.09 

13.1 

418.4 

23.85 

62.8 

451.7 

25.75 

94.7 

300.0 

42.54 

1120.07 

17.10 

444.0 

329.9 

18.80 

681.0 

257.0 

51.22 

1348.61 

14.65 

224.0 

283.5 

16.16 

371.0 

Side-on  Pressure  (kPa)  Side— on  Pressure  (k?a)  Sid 


Figure  14.  Hydrocode  Predictions  -  Short  Driver,  314  kPa 
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Impulse  (kPa— sae)  Impulse  (kPa— sec)  Impulse  OcPa-aec) 


Side— on  Pressure  (k Pa)  Side— on  Pressure  (kPa)  Side— on  Pressure  (kPa) 


Figure  16.  Hydrocode  Predictions  -  Short  Driver,  4413  kPa. 
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Impulse  (kPa-isc)  Impulse  (kPa-seo)  Impulse  (kPa-sec) 


Figure  17.  Hydroc.ode  Predictions  -  Short  Driver,  14480  kPa. 
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the  predicted  pressures  from  the  code  varied  both  above  and  below  the 
experimental  values.  The  general  decay  along  the  test  section  was  similar 
to  the  measured  experimental  values. 

The  positive  durations  predicted  were  5  to  10  percent  smaller  than  the 
experimental  records  showed.  Station  3  should  not  be  compared  here  (except 
for  peak  pressure)  since  it  was  near  the  open  end  of  the  computational 
shock  tube.  The  rarefaction  from  the  open  end  shortened  the  positive  dura¬ 
tion.  The  experiments  used  a  long  test  section  with  a  closed  end  so  no 
rarefactions  could  occur. 

The  agreement  in  negative  phase  between  the  code  and  experimental 
values  was  not  nearly  as  close.  The  code  values  of  negative  pressure  were 
anywhere  from  30  to  100%  lower  (Station  1,  7  diameters;  before  open  end 
rarefaction  arrived  at  station)  than  was  actually  measured.  Overall,  the 
code  predictions  were  quite  useful  to  predict  the  general  pressure-time 
waveforms  from  the  various  driver  configurations,  given  the  limitations 
discussed . 

3.3  Cold  Gas/Recompression  Fan  Effects.  Figures  18-22  show  some  of 
the  code  predictions  of  pressure-time  records  for  the  long  driver.  As  was 
noted  earlier,  side-on  overpressure  begins  to  be  predicted  too  high  for 
driver  pressures  of  about  3000  kPa.  Both  the  cold  gas  and  recompression 
fan  (backward  facing  shock)  arrival  are  very  apparent.  Table  5  summarizes 
and  compares  the  code  results  with  that  of  the  test  data.  Figures  23-25 
display  these  comparisons  in  graphical  form. 

Theoretically,  the  cold  gas  should  not  appear  on  side-on  overpressure 
records  at  all.  It  does,  however,  appear  on  the  experimental  records 
although  the  code  cannot  predict  it  for  the  side-on  records.  It  does,  of 
course,  predict  it  for  the  stagnation  records.  See  Figure  21.  Experi¬ 
mentally,  the  cold  gas  region  tends  to  diffuse  and  be  quite  turbulent  in 
nature.  The  code,  instead,  predicts  very  smooth  pressure  change  upon  the 
cold  gas  arrival;  and  it  predicts  a  slower  arrival  than  experimentally 
recorded.  The  cold  gas  average  drag  enhancement  is  predicted  by  the  code 
to  be  higher  than  the  data  suggests.  It  should  be  noted  that  higher  spikes 
do  exist  on  the  data  records.  Only  average  ratio  values  obtained  from  the 
impulse  data  were  used  to  plot  Figure  24  from  Table  5-C.  It  should  be 
noted  that  the  enhancement  is  about  constant  along  the  test  section  for  a 
given  driver  pressure. 

The  second  effect  that  changed  the  pressure-time  waveform  (and  hence 
yield)  is  shown  by  the  code  results  in  Figure  20.  Notice  the  very  sharp 
cutoff  of  the  positive  phase  and  the  relatively  low,  and  extensive  negative 
phase  of  the  records.  This  effect  shows  up  sooner  on  the  data  records  and 
is  more  of  a  rarefaction  decay  than  a  discontinuous  drop  as  the  code  sug¬ 
gests.  Experimentally,  the  cold  gas  and  recompression  fan  effects  move 
back  from  the  shock  front  with  travel  along  the  test  section.  The  data 
show  that  Station  70  (7  dia.)  cannot  be  used  with  a  cold  driver  operation 
above,  perhaps,  100  kPa  side-on  overpressure  at  the  station. 

A  shorter  driver,  and  shorter  duration  record,  ought  to  minimize  the 
cold  gas/recompression  fan  effects.  Figure  26  shows  a  comparison  of 
records  from  the  long,  288  cm,  driver  with  those  obtained  from  a  short, 
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Figure  18.  Hydrocode  Side-On  Overpressure  Predictions  - 
Long  Driver;  2034,  3103,  and  4220  kPa. 
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Side-on  Pressure  (kPa)  Sid 


Figure  18.  Hydrocode  Side-On  Overpressure  Predictions  - 

Long  Driver;  2034,  3103,  and  4220  kPa  (Cont'd). 
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Figure  19.  Hydrocode  Stagnation  Overpressure  Predictions  - 
Long  Driver;  2034,  3103,  and  4220  kPa. 
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B.  Station  at  20  Diameters 

Figure  19.  Hydrocode  Stagnation  Overpressure  Predictions  - 
Long  Driver;  2034,  3103,  and  4220  kPa  (Cont'd). 
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A.  Station  at  7  Diameters 


Figure  20.  Hydrocode  Side-On  Overpressure  Predictions  - 
Long  Driver;  7480,  10963,  and  13169  kPa. 
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B.  Station  at  20  Diameters 


Figure  20.  Hydrocode  Side-On  Overpressure  Predictions  - 

Long  Driver;  7480,  10963,  and  13169  kPa  (Cont'd) 
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Figure  21.  Hydrocode  Stagnation  Overpressure  Predictions  - 
Long  Driver;  7480,  10963,  and  13169  kPa. 
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Figure  21.  Hydrocode  Stagnation  Overpressure  Predictions  - 
Long  Driver;  7480,  10963,  and  13169  kPa  (Cont'd) 
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Figure  22.  Hydrocode  Dynamic  Overpressure  Predictions  - 
Long  Driver;  7480,  10963,  and  13169  kPa. 
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Station:  3 
<  0.40671  > 


Station:  3 
<  0.40671  > 


Figure  22.  Hydrocode  Dynamic  Overpressure  Predictions  - 

Long  Driver;  7480,  10963,  and  13169  kPa  (Cont'd) 
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TABLE  5.  COMPARISON  OF  COMPUTER  CODE  SIMULATION  WITH  MEASURED 
PARAMETERS;  288  CM  DRIVER,  16:1  THROAT  RATIO 


A.  Shock  Front  Overpressure,  kPa 


Shot 

Driver 

BRL-Q1D  Code 

Experiment 

Number 

Pressure 

Station 

Station 

kPa 

70 

150 

200  240 

.  280 

70 

150 

200 

240 

280 

18 

300 

13.1 

Same 

13.8 

— 

13.2 

13.3 

13.0 

17 

603 

24.4 

Same 

23.5 

22.8 

23.2 

23.0 

22.0 

19 

1138 

35.6 

Same 

38.0 

38.3 

38.0 

37.6 

37.0 

20 

2034 

57.5 

Same 

60.4 

60.0 

61.2 

60.4 

57.0 

16 

3103 

92.5 

Same 

88.0 

85.0 

86.5 

82.5 

82.5 

21 

4220 

120.0 

Same 

106.5 

100.0 

100.0 

100.0 

93.0 

14 

5171 

No  Run 

131.5 

122.5 

120.5 

111.8 

110.0 

15 

7480 

180.8 

Same 

155.0 

156.0 

153.8 

155.0 

150.0 

22 

10963 

238.5 

Same 

193.0 

193.0 

190.0 

191.0 

177.5 

23 

13169 

257.7 

Same 

215.0 

212.0 

212.0 

210.0 

203.0 

B •  Cold  Gas 

Arrival 

After  Shock 

Front 

,  ms 

20 

2034 

10.1 

— 

-  - 

— 

— 

— 

— 

— 

— 

16 

3103 

5.7 

— 

-  - 

— 

5.0 

— 

— 

— 

— 

21 

4220 

4.2 

9.2 

12.8 

— 

4.2 

— 

— 

— 

— 

14 

5171 

No  Run 

2.8 

7.3 

10.4 

— 

— 

15 

7480 

2.6 

6.3 

8.4  9.7 

11.6 

2.2 

6.1 

8.5 

10.7 

— 

22 

10963 

1.7 

4.0 

5.1  6.8 

8.0 

1.8 

4.7  • 

6.0 

7.6 

8.8 

23 

13169 

1.6 

3.7 

4.7  5.5 

6.6 

1.7 

4.3 

5.2 

6.6 

7.5 

TABLE  5.  COMPARISON  OF  COMPUTER  CODE  SIMULATION  WITH  MEASURED 
PARAMETERS;  288  CM  DRIVER,  16:1  THROAT  RATIO  (CONT'D) 


Shot  Driver 

Number  Pressure 


kPa 

20 

2034 

16 

3103 

21 

4220 

14 

5171 

15 

7480 

22 

10963 

23 

13169 

21 

4220 

14 

5171 

15 

7480 

22 

10963 

23 

13169 

C.  Cold  Gas  Drag  Enhancement  Ratio 
BRL-Q1D  Code 
Station 


Experiment 

Station 


70 

150 

200 

240 

280 

70 

150 

200 

240 

1.33 

— 

— 

— 

— 

— 

— 

— 

1.51 

— 

— 

— 

— 

— 

— 

1.75 

1.77 

— 

— 

— 

— 

— 

— — 

No  Run 

— 

1.32 

1.20 

1.15 

2.42 

2.27 

2.27 

2.22 

— 

— 

1.48 

1.50 

1.50 

3.04 

2.95 

2.95 

2.90 

2.92 

— 

1.58 

1.60 

1.60 

3.33 

3.48 

3.45 

3.45 

3.36 

— 

1.80 

1.90 

1.90 

D.  Recompression 

Fan  Arrival  After 

Shock 

Front,  ms 

_ 

— 

— 

— 

— 

4.2 

— 

— 

— — 

No  Run 

3.5 

— 

— 

— 

8.0 

17.6 

29.5 

— 

— 

3.5 

14.0 

— 

— 

5.3 

11.3 

15.5 

20.0 

25.2 

3.2 

10.0 

— 

— 

4.5 

9.8 

12.6 

16.0 

19.7 

3.8 

10.1 

12.9 

18.0 

COLD  GAS  DRAG  ENHANCEMENT  RATIO 


Figure  24.  Drag  Enhancement  Ratio  -  Caused  by  Cold 
Gas  Arrival  at  Test  Station. 
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SINGLE  DRIVER  -  10.16cm  DIA.*288cm  LONG  , 
16N  TEST  SECTION  TO  THROAT  AREA  RATIO 


TEST:  LBTS/57 


TIME,  MSEC 


TEST:  LBTS/57 


TIME,  MSEC 


TEST:  LBTS/57 


TIME,  MSEC 

C.  Short  Driver  -  Station  70,  16:1  Throat  Ratio 


Figure  26.  Comparison  of  Cold  Gas  Effects  for 

Different  Driver  Configurations  (Cont'd). 
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.  KPA-MSEC  IMPULSE.  KPA-KSEC  IMPULSE.  KPA-KSEC 


TEST:  LBTS/57 


TEST:  LBTS/57 


TEST:  LBTS/57 


D.  Short  Driver  -  Station  200,' 16:1  Throat  Ratio 


Figure  26.  Comparison  of  Cold  Gas  Effects  for  Different 
Driver  Configurations  (Cont'd). 
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.  KPA-MSEC  IMPULSE.  KPA-MSEC  IMPULSE.  KPA-MSEC 


TIME,  MSEC 


F.  Short  Driver  -  Station  200,  S3 

Figure  26.  Comparison  of  Cold  Gas  Off 
Driver  Configurations  (Con 


11.13  cm,  driver.  Indeed,  the  shorter  driver  does  show  the  lesser  effects 
on  the  records  but  are  still  unacceptable  for  a  correct  simulation  of  field 
yields  at  this  pressure  level  if  taken  at  7  diameters.  Station  200  (20 
dia.)  does  show  an  improvement  with  acceptable  records  as  seen  in  Figure 
26-D  and  26-F. 

Accordingly,  the  full-size  predictions  for  a  single,  cold  gas  driver 
simulator  as  presented  in  the  next  section  are  generally  for  the  20 
diameter  station. 

4.  ANALYSIS 

The  analysis  includes  a  discussion  of  the  full-size  simulator  yield 
predictions  as  calculated  from  side-on  overpressure  impulse  scaled  from 
Reference  8  data,  and  the  effects  of  changing  the  throat  baffle  ratio  at 
the  diaphragm  section. 

4.1  Yield  Predictions  as  a  Function  of  Driver  Configurations.  As 
shown  in  the  Test  Matrix  of  Table  1,  a  series  of  drivers  from  6.05  cm  to 
288  cm  was  used  on  the  model  for  a  driver  pressure  range  of  about  300  kPa 
to  1800  kPa.  Side-on  overpressures  at  7  diameters  varied  from  about  10  kPa 
to  225  kPa . 

Reference  9  presents  scaling  laws  relating  the  blast  parameters  between 
explosive  yields.  The  treatment  of  yield  here  will  consider  ground  bursts 
and  the  impulse  available  at  a  given  overpressure  from  a  certain  yield  in 
kilotons  (kt). 

The  calculations  of  predicted  yield  for  the  full-size  simulator  were 
made  by  first  scaling  up  the  impulse  from  the  model  experiments,  by  multi¬ 
plying  by  the  scale  factor  of  57.  At  the  same  time,  driver  configurations 
were  scaled  in  the  same  way.  Table  6  lists  the  full-size  predicted  peak 
overpressure  and  impulse  as  a  function  of  driver  pressure  and  configura¬ 
tion.  Predicted  records  of  pressure  versus  time  and  impulse  versus  time 
for  stations  at  7  to  20  diameters  for  the  full-size  simulator  are  shown  in 
Figures  27  and  28.  Figures  29  and  30  summarize  the  full-scale  impulse  as  a 
function  of  driver  parameters.  Driver  lengths  and  volumes  are  noted  on 
each  curve. 

Values  of  full-scale  impulse  from  Table  6  were  then  compared  to  the 
values  given  in  Figure  31.  This  figure  is  a  plot  of  data  taken  from  the 
Defense  Nuclear  Agency  (DNA)  blast  program  listed  in  Reference  8.  The 
yields  were  then  found  from  Equation  1. 

VT2  =  (Vw2)i/3>  (1) 

where  1^  is  the  impulse  for  a  given  reference  yield,  W^,  (in  this  case, 

1  kt)  and  I^  is  the  impulse  at  the  same  pressure  level  for  the  predicted 
yield,  W^.  Impulse  is  given  in  kPa-s  and  yield  in  kt  for  simulation  of 
nuclear  explosions  . 

It  should  be  noted  that  for  other  field  scenerios,  other  yields  may  be 
calculated  for  height-of-bursts  (HOB),  or  for  dynamic  overpressure  impulse, 
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TABLE  6.  PREDICTED  YIELD  FOR  FULL-SIZE  FROM  EXPERIMENTAL 
SIDE-ON  OVERPRESSURE  IMPULSE 


Shot 

Driver 

Driver 

Driver 

Throat 

Side-on 

Number- 

Pressure 

Length 

Volume 

Baffle  Ratio 

Overpressure 

Impulse 

Yield 

Station 

kPa 

m 

m3 

kPa 

kPa-ms 

kt 

4-200 

314 

6.34 

166.9 

16:1 

17.0 

1.30 

0.1 

47-200 

317 

19.37 

510.0 

18.6 

3.20 

1.9 

53-200 

338 

51.22 

1348.6 

14.0 

7.47 

34.0 

55-200 

359 

101.93 

2683.8 

13.0 

15.55 

365.0 

57-200 

352 

135.19 

3559.5 

13.0 

18.66 

632.0 

6-200 

1124 

6.34 

166.9 

34.1 

4.41 

1.3 

49-200 

1096 

19.37 

510.0 

49.1 

9.62 

7.6 

79-200 

1124 

41.09 

1081.9 

42.5 

20.50 

91.0 

78-200 

1117 

62.81 

1653.8 

40.5 

28.84 

178.0 

80-200 

1145 

83.07 

2187.2 

45.0 

41.04 

677.0 

8-200 

3137 

6.34 

166.9 

66.7 

10.57 

6.3 

69-200 

3137 

19.37 

510.0 

94.6 

24.02 

49.0 

70-200 

5171 

19.37 

510.0 

127.0 

36.54 

132.0 

75-200 

3137 

41.09 

1081.9 

94.0 

48.88 

411.0 

77-200 

3123 

62.81 

1653.8 

99.0 

72.66 

1275.0 

11-200 

8274 

6.34 

166.9 

113.4 

23.95 

44.0 

71-200 

8486 

19.37 

510.0 

172.4 

56.68 

419.0 

72-200 

11549 

19.37 

510.0 

202.5 

75.09 

882.0 

13-200 

14479 

6.34 

166.9 

155.8 

40.52 

67.0 

66-200 

13962 

10.68 

281.2 

192.0 

60.00 

466.0 

65-70 

17582 

3.44 

90.6 

200.0 

25.98 

37.0 

65-200 

17582 

3.44 

90.6 

141.0 

30.78 

74.0 

TABLE  6.  PREDICTED  YIELD  FOR  FULL-SIZE  FROM  EXPERIMENTAL 
SIDE-ON  OVERPRESSURE  IMPULSE  (CONT.) 


Shot 

Driver 

Driver 

Driver 

Number- 

Pressure 

Length 

Volume 

Station 

kPa 

m 

m3 

81-200 

483 

19.37 

510.0 

85-200 

483 

38.19 

1005.5 

89-200 

496 

54.12 

1425.0 

91-200 

545 

83.07 

2187.2 

35-200 

2027 

6.34 

166.9 

82-200 

1827 

19.37 

510.0 

85-200 

1827 

38.19 

1005.5 

90-200 

1834 

54.12 

1425.0 

37-200 

5192 

6.34 

166.9 

83-200 

5199 

19.37 

510.0 

87-200 

5240 

38.19 

1005.5 

40-70 

14479 

6.34 

166.9 

40-200 

14479 

6.34 

166.9 

84-200 

14789 

19.37 

510.0 

94-200 

903 

42.54 

1120.1 

Throat 
Baffle  Ratio 

33:1 


Side-on 

Overpressure 

Impulse 

Yield 

kPa 

kPa-ms 

kt 

15.7 

4.81 

7.3 

15.7 

9.09 

51.0 

13.5 

11.51 

138.0 

15.0 

20.69 

641.0 

39.0 

6.84 

4.0 

39.3 

16.74 

58.0 

40.6 

31.35 

356.0 

39.8 

41.04 

771.0 

83.4 

15.45 

14.9 

88.9 

37.60 

198.0 

94.9 

74.10 

1415.0 

195.3 

32.37 

72.0 

142.4 

36.28 

121.0 

200.0 

86.96 

1386.0 

13.0 

18.80 

681.0 

64:1 


PRESSURE,  KPA  PRESSURE,  KPA 
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IMPULSE.  KPA-MSEC  IMPULSE.  KPA-MSEC  IMPULSE.  KPA-MSEC  3WJLSE.  KPA-MSEC 


PRESSURE,  KPA 
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KPA-HSEC  IMPULSE.  KPA-MSEC  IMPULSE.  KPA-MSEC 


!.  KPA-HSEC  IMPULSE.  KPA-NSEC  IMPULSE.  KPA-MSEC 


FULL-SCALE  OVERPRESSURE  IMPULSE  ,  k  P 


SINGLE  DRIVER,  5.79m  DIA. 

33:1  TEST  SECTION  TO  THROAT  AREA  RATIO 


Figure  30.  Full-Scale  Impulse  vs  Side-On  Overpressure  at  20  Diameters; 
for  Different  Driver  Lengths,  33:1  Throat  Ratio. 
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if  desired.  The  method  of  yield  calculation  would  be  the  same,  but 
impulses  would  be  compared  to  another  impulse-yield  curve  instead  of  Figure 
31. 


The  yields  as  a  function  of  side-on  overpressure  and  driver/ throat 
ratio  configuration  are  presented  in  Figures  32  and  33.  Data  are.  plotted 
from  both  Tables  3  and  6  in  order  to  compare  Stations  7  and  20.  These 
curves  are  ordered  very  well  except  for  points  on  the  223  kPa  set  at  the 
high  yield  end.  There,  the  recompression  fan  has  shortened  the  positive 
wave  duration  and  the  corresponding  yield  for  a  particular  driver  con¬ 
figuration.  As  a  result,  the  curve  turns  away  from  the  left  axis  to  the 
right  -  indicating  a  longer  driver  is  needed. 

Table  6,  above,  lists  the  calculated  full-size  yields  in  the  last 
column,  as  a  function  of  the  various  test  station  pressure  levels  and  the 
driver  configurations  used  in  the  tests.  The  table  illustrates  the  range 
of  yields  available  for  a  single,  cold  gas  driver  design.  Most  are 
adequate  except  at  low  yields  (below  50  kt)  for  test  station  pressures 
above  200  kPa.  The  decay  rate  caused  by  the  relatively  short  positive 
durations  needed,  will  be  substantial.  The  driver  pressure  then  would  have 
to  be  increased  to  achieve  the  desired  test  station  pressure  level.  This 
would  add  additional  expensive  material  design  constraints  on  the  driver. 

A  heated  driver  design  would  probably  be  a  better  trade-off  for  a  large 
simulator  than  designing  for  the  expected  increased  driver  pressure  needed. 

4.2  Throat  Baffle  Effects.  A  major  element  in  the  design  of  the 
driver  configuration  is  the  test  section  to  throat  baffle  area  ratio 
(baffle  ratio).  As  noted  in  the  tables  and  graphs,  throat  ratios  of  16:1, 
33:1,  and  64:1  were  listed  for  this  set  of  experiments.  The  throat  baffle 
size  is  most  important  when  designing  for  large  yields  (600  kt/at  low  test 
pressures  (13-15  kPa  in  these  tests)).  Values  from  Table  6  show  the  effect 
of  the  change  in  throat  baffle  ratio.  Figure  34  illustrates  the  similarity 
of  waveforms  obtained  by  varying  the  driver  length,  the  driver  pressure, 
and  the  baffles  ratio.  To  summarize,  the  increasing  baffle  ratio  (smaller 
area  baffle)  restricts  the  driver  flow  to  add  more  duration  to  a  waveform 
produced  by  a  given  driver  length.  The  peak  overpressure  produced  becomes 
a  smaller  value  at  the  same  time. 

Figure  34  illustrates  the  results  of  three  baffle  ratios  used.  The 
maximum  driver  length  needed  at  the  13-15  kPa  overpressure  was  42.5  m 
(full-size)  when  it  was  used  with  the  64:1  baffle  ratio;  whereas,  a  135 
meter  length  is  needed  for  the  16:1  baffle  ratio.  A  penalty  in  increased 
driver  pressure  must  be  paid,  however,  for  the  shortened  driver  as  is  seen 
from  Figure  35.  The  driver  pressure  (noted  in  parenthesis)  must  be  roughly 
doubled  for  a  throat  ratio  increase  of  two  times  if  the  test  section  pres¬ 
sure  is  to  be  maintained.  At  low  driver  pressures,  the  smaller  length 
driver  probably  will  still  be  a  cost  effective  design  for  the  full-size 
simulator . 

5.  SUMMARY  AND  CONCLUSIONS 

Experimental  data  have  been  obtained  and  compared  with  BRL-Q1D  hydro¬ 
code  predictions  of  blast  effects  obtained  for  a  1:57  scale  shock  tube 
model  of  a  single  cold  air  driver  for  large  blast  thermal  simulator. 
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The  text  -is  continued  on  page  84. 


PRESSURE,  KPA  PRESSURE,  KPA 


TEST  SECTION  TO  THROAT  BAFFLE  AREA  RATIO 


SINGLE  DRIVER 
STATION  20 DIA. ,  640  kT 


DRIVER  LENGTH,  meters 


Figure  35.  Effect  of  Throat  Baffle  on  Driver  Length. 
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Thermal  effects  were  not  a  part  of  this  study.  BRL  has  the  task  of  prepar¬ 
ing  a  design  for  the  Army  of  a  large  simulator  capable  of  full-size  multi¬ 
service  equipment  tests  to  simulate  the  blast  and  thermal  effects  from 
nuclear  events.  The  experiments  and  code  comparisons  reported  here  were 
designed  to  give  insight  into  a  number  of  blast  parameters  fundamental  to 
the  design  of  such  a  simulator. 

Parameters  studied  were:  blast  wave  decay,  overpressure  in  the  test 
section  as  a  function  of  driver  pressure,  blast  wave  duration  and  impulse 
as  a  function  of  driver/throat  baffle  configuration,  and  effects  of  the 
cold  air  driver  gas/recompression  fan  arrival  at  a  particular  test  station. 

It  was  concluded  that  a  station  at  20  test  section  diameters  was  needed 
to  minimize  the  cold  gas/recompression  fan  effects,  at  test  station  pres¬ 
sures  of,  and  above,  100  kPa.  A  pressure  decay  problem  was  also  evident 
above  this  pressure  for  low  yields  (short  duration,  blast  waves).  This 
could  be  corrected  with  the  use  of  still  higher  driver  pressures.  The  sug¬ 
gested  preliminary  large  simulator  design  had  not  provided  for  these  higher 
pressures  . 

A  heated  driver  to  replace  the  cold  driver  design  should  probably  be 
considered.  Nearer  test  stations  might  possibly  be  utilized  in  the  heater 
driver  design.  Acceptable  levels  of  recompression  fan  and  peak  pressure 
decay  effects,  hopefully,  could  be  obtained  with  such  a  design.  Cold  gas 
effects  would,  of  course,  have  been  eliminated  by  the  heated  driver. 

Throat  baffles  would  still  be  needed  to  minimize  the  driver  length  at  the 
low  pressure/high  yield  end  of  the  simulator's  operating  range. 

The  BRL-Q1D  code  was  quite  helpful  in  predicting  the  general  waveforms 
for  the  various  driver/baffle  configurations  tested.  The  predicted 
enhancement  of  the  drag  because  of  cold  gas/recompression  fan  effects  from 
the  code  were  higher  than  actually  measured.  Also,  the  code  predicted  more 
test  time  to  be  available  before  the  cold  gas  or  recompression  fan  arrived 
at  a  given  test  station.  The  Q1D  code  did  furnish  some  helpful  guidelines 
for  a  very  complicated  blast  parameter  study;  it  did  it  efficiently  and  at 
a  low  cost.  Further  computations  utilizing  two-dimensional  hydrocodes, 
accounting  for  the  real  gas  effects  present  experimentally,  would  of 
course,  provide  better  experimental  agreement.  This  would  be  done  at 
greater  costs  in  both  computing  time  and  effort,  however. 
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Lawrence  Livermore  Laboratory 
Technical  Information  Division 
P.  0.  Box  808 
Livermore,  CA  94550 

1  Director 

Los  Alamos  Scientific  Lab 
ATTN:  Dr.  J.  Taylor 
P.  0.  Box  1663 
Los  Alamos,  NM  87544 

2  Director 

Sandia  National  Laboratories 
ATTN:  Info  Dist  Div 

Dr.  W.  A.  von  Riesemann 
(Div  6442) 

Albuquerque,  NM  87115 
1  Director 

National  Aeronautics  and  Space 
Administration 
George  C.  Marshall  Space 
Flight  Center 
Huntsville,  AL  35812 

Director 

National  Aeronautics  and  Space 
Administrtion 

Aerospace  Safety  Research  and 
Data  Institute 
21000  Brook  Park  Road 
Lewis  Research  Center 
Cleveland,  OH  44135 
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DISTRIBUTION  LIST 


No .  of 

Copies  Organization 

i  Director 

National  Aeronautics  and  Space 
Administration 
Scientific  and  Technical 
Information  Facility 
P.  0.  Box  8757 
Baltimore /Washington 
International  Airport, 

MD  21240 

1  National  Academy  of  Science 
ATTN:  Mr.  D.  G.  Groves 
2101  Constitution  Avenue,  NW 
Washington,  DC  20418 

10  Central  Intelligence  Agency 
OIR/DB/Standard 
GE47  HO 

Washington,  DC  20505 

1  DNA  Information  and  Analysis 
Center 
Kaman  Tempo 
ATTN:  DA SI AC 

816  State  Street 
P.0.  Drawer  00 
Santa  Barbara,  CA  93102 

1  Institute  of  Makers  of 
Explosives 

ATTN:  Mr.  F.  P.  Smith,  Jr., 
Executive  Director 
1575  Eve  St.,  N.W.,  Suite  550 
Washington,  DC  20005 

1  Aberdeen  Research  Center 
ATTN:  Mr.  John  Keefer 
30  Diamond  St. 

P.  0.  Box  548 
Aberdeen,  MD  21001 

1  Agbabian  Associates 

ATTN:  Dr.  D.  P.  Reddy 

250  N.  Nash  Street 
El  Segundo,  CA  90245 


No .  of 

Copies  Organization 

1  Ammann  &  Whitney 
ATTN:  Mr.  N.  Dobbs 
Suit  1700 

Two  World  Trade  Center 
New  York,  NY  10048 

1  Black  &  Veatch  Consulting 
Engineers 

ATTN:  Mr.  H.  L.  Callahan 
1500  Meadow  Lake  Parkway 
Kansas  City,  MO  64114 

1  Dr.  Wilfred  E.  Baker 

Wilfred  Baker  Engineering 

P.  0.  Box  6477 

San  Antonio,  TX  78209 

1  Aeronautical  Research 

Associates  of  Princeton, 
Inc . 

ATTN:  Dr.  C.  Donaldson 
50  Washington  Road, 

P.O.Box  2229 
Princeton,  NJ  08540 

1  Applied  Research  Associates, 
Inc . 

ATTN:  Mr.  J.  L.  Drake 
1204  Openwood  Street 
Vicksburg,  MS  39180 

1  J.  G.  Engineering  Research 

Associates 
3831  Menlo  Drive 
Baltimore,  MD  21215 

2  Kaman-AviDyne 

ATTN:  Dr.  N.P. Hobbs 

Mr.  S.  Criscione 
Northwest  Industrial  Park 
83  Second  Avenue 
Burlington,  MA  01803 

3  Kaman-Nuclear 

ATTN:  Dr.  F.  H.  Shelton 
Dr.  D.  Sachs 
Dr.  R.  Keffe 

1500  Garden  of  the  Gods  Road 
Colorado  Springs,  CO 
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1  Knolls  Atomic  Power  Laboratory 
ATTN:  Dr.  R.  A.  Powell 
Schenectady,  NY  12309 

.1  McDonnell  Douglas  Astronautics 
Western  Division 
ATTN:  Dr.  Lea  Cohen 
5301  Bosla  Avenue 
Huntington  Beach,  CA  92647 

1  Physics  International 
2700  Merced  Street 
San  Leandro,  CA  94577 

1  R&D  Associates 

ATTN:  Mr.  John  Lewis 

P.  0.  Box  9695 

Marina  del  Rey,  CA  90291 

1  R&D  Associates 
ATTN:  G.  P.  Ganong 
P.  0.  Box  9335 
Albuquerque,  NM  87119 

2  The  Boeing  Company 
Aerospace  Division 
ATTN:  Dr.  Peter  Grafton 

Dr.  D.  Strome 
Mail  Stop  8C-68 
P.  0.  Box  3707 
Seattle,  WA  98124 

2  AVCO  Corporation 

Structures  and  Mechanics  Dept. 
ATTN:  Dr.  William  Eroding 
Dr.  J.  Gilmore 
201  Lowell  Street 
Wilmington,  MA  01887 

1  Aerospace  Corporation 
P.0.  Box  92957 
Los  Angeles,  CA  90009 

1  General  American 

Transportation  Corp. 

General  American  Research  Div. 
ATTN:  Dr.  J.  C.  Shang 
7449  N.  Natchez  Avenue 
Niles,  IL  60648 


No .  of 

Copies  Organization 

1  Hercules,  Inc. 

ATTN:  Billings  Brown 
Box  93 

Magna,  UT  84044 

1  Mason  &  Hanger-Silas  Mason 
Co.,  Inc. 

Plantex  Plant 
P.  0.  Box  647 
Amarillo,  TX  79117 

1  Lovelace  Research  Institute 
ATTN:  Dr.  E.  R.  Fletcher 
P.  0.  Box  5890 
Albuquerque,  NM  87115 

1  Massachusetts  Institute  of 
Technology 

Aeroelastic  and  Structures 
Research  Laboratory 
ATTN:  Dr.  E.  A.  Witmar 
Cambridge,  MA  02139 

1  Monsanto  Research  Corporation 
Mound  Laboratory 
ATTN:  Frank  Neff 
Miamisburg,  OH  45342 

1  Science  Applications,  Inc. 

Suite  310 

1216  Jefferson  Davis  Highway 
Arlington,  VA  22202 

2  Battelle  Memorial  Institute 
ATTN:  Dr.  L.  E.  Hulbert 

Mr.  J.  E.  Backofen,  Jr. 
505  King  Avenue 
Columbus,  OH  43201 

1  Georgia  Institute  of  Tech 
ATTN:  Dr.  S.  Atluri 
225  North  Avenue,  NW 
Atlanta,  GA  30332 

1  IIT  Research  Institute 
ATTN:  Mrs.  H.  Napadensky 
10  West  35  Street 
Chicago,  IL  60616 
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2  Southwest  Research  Institute 
ATTN:  Dr.  H.  N.  Abramson 
Dr.  U.  S.  Lindholm 
8500  Culebra  Road 
San  Antonio,  TX  78228 

1  Brown  University 

Division  of  Engineering 
ATTN:  Prof.  R.  Clifton 
Providence,  RI  02912 

1  Florida  Atlantic  University 
Dept.,  of  Ocean  Engineering 
ATTN:  Prof.  K.  K.  Stevens 
Boca  Raton,  FL  33432 

1  Texas  A&M  University 
Department  of  Aerospace 
Engineering 

ATTN:  Dr.  James  A.  Stricklin 
College  Station,  TX  77843 

1  University  of  Alabama 
ATTN:  Dr.  T.  L.  Cost 
P.  0.  Box  2908 
University,  AL  35486 


No .  of 

Copies  Organization 

1  University  of  Delaware 

Department  of  Mechanical  and 
Aerospace  Engineering 
ATTN:  Prof.  J.  R.  Vinson 
Newark,  DE  19711 


Aberdeen  Proving  Ground 

Dir,  USAMSAA 
ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 
Cdr ,  USATECOM 

ATTN:  AMSTE-SI-F 
Cdr,  CRDC,  AMCC0M 
ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-SPS-IL 
Cdr,  US  Army  Toxic  and 

Hazardous  Materials  Agency 
ATTN:  AMXTH-TE 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Your  comments/answers  to  the  items/questions  below  will 
aid  us  in  our  efforts. 

1 .  BRL  Report  Number _ Date  of  Report 

2.  Date  Report  Received _ 


3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or 
other  area  of  interest  for  which  the  report  will  be  used.) 


4.  How  specifically,  is  the  report  being  used?  (Information  source,  design 
data,  procedure,  source  of  ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far 
as  man-hours  or  dollars  saved,  operating  costs  avoided  or  efficiencies  achieved, 
etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future 
reports?  (Indicate  changes  to  organization,  technical  content,  format,  etc.) 


Name 


CURRENT 

ADDRESS 


Organization 


Address 


City,  State,  Zip 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the 
New  or  Correct  Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


OLD 

ADDRESS 


Name 


Organization 


Address 


City,  State,  Zip 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


—  -  -  -  -  - - -  FOLD  HERE  - 

Director 

US  Army  Ballistic  Research  Laboratory 
ATTN:  DRXBR-OD-ST 

Aberdeen  Proving  Ground,  MD  21005-5066 


NO  POSTAGE 
NECESSARY 
IF  MAILED 
IN  THE 

UNITED  STATES 


OFFICIAL  BUSINESS 

PENALTY  FOR  PRIVATE  USE,  *300 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  NO  12062  WASHINGTON, DC 
POSTAGE  WILL  BE  PAIO  BY  DEPARTMENT  OF  THE  ARMY 


Director 

US  Army  Ballistic  Research  Laboratory 
ATTN:  DRXBR-OD-ST 

Aberdeen  Proving  Ground,  MD  2100S-9989 


FOLD  HERE 


